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Phase separation is a ubiquitous process in seafloor hydrothermal vents, creating a large 
range of salinities. Toxic elements (e.g., arsenic) partition into the vapor phase, and thus 
can be enriched in both high and low salinity fluids. However, investigations of microbial 
diversity at sites associated with phase separation are rare. We evaluated prokaryotic 
diversity in arsenic-rich shallow-sea vents off Milos Island (Greece) by comparative analysis 
of 16S rRNA clone sequences from two vent sites with similar pH and temperature 
but marked differences in salinity. Clone sequences were also obtained for aioA-Wke 
functional genes (AFGs). Bacteria in the surface sediments (0-1.5 cm) at the high salinity 
site consisted of mainly Epsilonproteobacteria (Arcobacter sp.), which transitioned to 
almost exclusively Firmicutes (Bacillus sp.) at ~10cm depth. However, the low salinity 
site consisted of Bacteroidetes (Flavobacteria) in the surface and Epsilonproteobacteria 
(Arcobacter sp.) at ~10cm depth. Archaea in the high salinity surface sediments 
were dominated by the orders Archaeoglobales and Thermococcales, transitioning to 
Thermoproteales and Desulfurococcales (Staphylothermus sp.) in the deeper sediments. 
In contrast, the low salinity site was dominated by Thermoplasmatales in the surface and 
Thermoproteales at depth. Similarities in gas and redox chemistry suggest that salinity 
and/or arsenic concentrations may select for microbial communities that can tolerate 
these parameters. Many of the archaeal 16S rRNA sequences contained inserts, possibly 
introns, including members of the Euryarchaeota. Clones containing AFGs affiliated with 
either Alpha- or Betaproteobacteria, although most were only distantly related to published 
representatives. Most clones (89%) originated from the deeper layer of the low salinity, 
highest arsenic site. This is the only sample with overlap in 16S rRNA data, suggesting 
arsenotrophy as an important metabolism in similar environments. 
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INTRODUCTION 

Geochemical gradients often dictate microbial community struc- 
tures, and metabolic processes directly influence these gradients. 
In hydrothermal systems, steep gradients generate redox disequi- 
libria that can provide the necessary energy for a wide diversity 
of Archaea and Bacteria. Since their discovery in the 1970s, deep- 
sea hydrothermal environments have garnered much attention in 
this regard, with a number of studies linking geochemical com- 
position, reaction energetics, and metabolic diversity (e.g., Flores 
et al., 2011; Meyer-Dombard et al., 2011). Their shallow-sea 
counterparts provide highly complementary research opportu- 
nities that are often overlooked. These systems are ubiquitous, 
readily accessible, and geochemically diverse; many exhibit geo- 
logical, chemical, and biological characteristics similar to those 
found at deep-sea vents. 



Shallow-sea and deep-sea hydrothermal systems also differ in 
several key ways (Dando et al., 1999; Pichler, 2005; Tarasov et al, 
2005; Price et al., 2012). For example, the evolved discharging 
hydrothermal fluids in shallow-sea systems may have originally 
been seawater, meteoric water, or a mixture of the two. They 
often occur within the photic zone and thus provide the oppor- 
tunity for both chemosynthetic and photosynthetic microbial 
metabolisms. They are most often associated with arc volcan- 
ism, which provides the heat source that drives hydrothermal 
circulation. Well known examples of shallow-sea hydrothermal 
systems are found near the Tabar-Feni (Pichler et al., 1999, 
2006), the Aeolian (Italiano and Nuccio, 1991; Amend et al, 
2003), the South Aegean (or Hellenic; Varnavas and Cronan, 
1988; Dando et al, 1995), the Caribbean (McCarthy et al., 2005), 
and the Kurile-Kamchatka island arcs (Tarasov et al, 1990). 
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FIGURE 1 | (A) Location of Milos Island (box) and other calc-alkaline 
volcanoes (shaded) along the Aegean island arc (dotted line). (B) Milos 
Island and the location of Palaeochori and Spathi Bays (box). Stippled 
offshore areas around the island are mapped gas emissions by echo 
sounding (Dando et al., 1995). (C) Location of our study site (box) in 
Palaeochori Bay. 



In addition, some shallow-sea vent systems also occur in other 
areas related to faulting (e.g., Bahia Conception; Forrest et al., 
2005) or serpentinization reactions (e.g., New Caledonia; Cox 
et al., 1982). Many of these shallow-sea systems are sediment 
covered, allowing steep vertical and horizontal geochemical gra- 
dients to evolve as reduced, low pH, high temperature fluids 
mix with overlying seawater (e.g., Price et al, 2007). Across 
the geochemical gradients, "niches" of potential energy develop 
that constrain the resident microbial communities. While most 
shallow-sea hydrothermal fluids seem to be sulfur- and/or iron- 
rich, they can also be elevated in potentially toxic elements 
such as As, Sb, Cr, Pb, Cd, Cu, Zn, and Hg (Varnavas and 
Cronan, 1988; Koski et al., 2001; Pichler et al., 2006; Price et al., 
2012). 

Many seafloor hydrothermal systems are modified by phase 
separation (German and Von Damm, 2003; Von Damm et al, 
2003), resulting in vent fluids that can vary drastically in salin- 
ity, from <6 to ~200% of seawater values (German and Von 
Damm, 2003). The resultant high and low salinity fluids differ 
markedly in the concentrations of many solutes, but the rela- 
tion to microbial community structures has not been adequately 
considered. Previous investigations indicate that the vapor phase 
fluids maybe enriched in important electron donors, such as H2, 
relative to the high salinity brine phase, and may contribute to 
stratification of microbial communities (Nakagawa et al., 2005; 
Nunoura and Takai, 2009). Arsenic is known to partition into 
the vapor phase (Pokrovski et al., 2002; Price et al., 2012). If 
the hydrothermal reservoir is enriched in arsenic and undergoing 
phase separation, this element can therefore be elevated in both 
high and low salinity discharging hydrothermal fluids. The high- 
est As levels in a marine hydrothermal fluid were reported for the 
shallow-sea hydrothermal vent system in Palaeochori Bay, Milos 
Island, Greece (Price et al, 2012). There, high arsenic concen- 
trations were reported for high and low salinity fluids; typically 
in the range of 30 u,M, but as high as 78 \iM in the low salin- 
ity fluids, suggesting As partitions into the vapor phase in this 
system. 

The purpose of this investigation was to identify the dominant 
bacterial and archaeal lineages in Palaeochori Bay in an attempt 
to understand the link between geochemical gradients and micro- 
bial groups, including arsenic oxidizers. These relationships were 
investigated at two sites with similar pH and temperature, but 
different salinities. 

SITE CHARACTERISTICS 

Gas discharge defines an ~35km 2 area of hydrothermal venting 
around Milos, making it one of the largest shallow-sea hydrother- 
mal systems described to date (Figures 1A,B; Dando et al., 1995). 
The most intense venting occurs in Palaeochori Bay, where 
CC^-rich gases and hydrothermal fluids discharge through sand 
(Figures 1A-C, 2). Nearby — but little explored — Spathi Bay, to 
the east, also features abundant hydrothermal venting, although 
free gas emissions are less abundant. The high salinity hydrother- 
mal fluids in Palaeochori Bay were slightly acidic (pH ~5), 
sometimes highly sulfidic (up to 3 mM), and warm (ambient to 
~110°C; Stuben and Glasby, 1999; Valsami- Jones et al., 2005; 
Price et al, 2012). 



The occurrence of orange, white, or brown microbial mats and 
hydrothermal precipitates on the seafloor sediments are a major 
feature of the hydrothermal system. The different colors corre- 
late to different sediment temperatures, where the hottest areas 
(>90°C) are typically covered by bright yellow native sulfur or 
orange-colored arsenic sulfide deposits (Figure 2A). White mats 
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FIGURE 2 | (A) Photograph of the RP site investigated in this study. P-0, 
Palaeochori orange; P-W, Palaeochori white; and P-B, Palaeochori 
background. Geochemical profiles and sediment cores were sampled at 
~P-W. (B) Photograph of the TW site. Geochemical and microbiological 
samples for this study are indicated at the location of TW. 



make up the largest areas of hydrothermal venting throughout the 
bay, and display mid-range temperatures (e.g., ~45-85°C). In the 
lower temperature areas (~30-35°C), brown colored manganese 
and/or iron oxide deposits are present (Figure 2A). 

Several microbiological investigations were carried out for 
the high salinity sites in Palaeochori Bay, although none to 
date targeted low salinity environments. They included culture- 
dependent and molecular-based microbiological studies, as well 
as the isolation of several new bacterial and archaeal species 
(Jochimsen et al., 1997; Dando et al., 1998; Brinkhoff et al., 
1999; Sievert et al, 1999, 2000a,b; Arab et al, 2000; Sievert 
and Kuever, 2000; Schlesner et al., 2001; Bayraktarov et al, 
2013). One study investigated the microbial ecology and arsenic 
functional genes in lower temperature 0-6 cm pooled sedi- 
ments at the same high salinity site investigated here (Nitzsche, 
2010). Results indicated a high diversity in the bacterial com- 
munity sequences of Gamma-, Delta-, and Epsilonproteobacteria, 
Bacteroidetes, and Cyanobacteria. Closest cultivated relatives were 
associated with photo- and chemotrophic lifestyle and seemed 



mostly to be involved in sulfur cycling. Archaeal sequences pre- 
dominantly affiliated with Thermoplasmatales or Thermococcales. 
Screening for arsenic functional genes indicated affiliation with 
either Alpha- or Betaproteobacteria, although affiliation with 
Gammaproteobacteria was noted. These studies indicated that 
metabolisms based on sulfur, either oxidation of H2S or reduc- 
tion of SO^, may be dominant at the high salinity sites, 
although dissimilatory iron reducing (DIR) Bacteria are abun- 
dant in some Felll-rich sediments. They also suggest that both 
heterotrophic and chemotrophic metabolisms are possible, with 
a range from mesophilic to hyperthermophilic functioning, and 
that arsenite oxidase genes are present at the high-salinity, lower 
temperature sites. This study expands on these investigations 
by studying the microbial ecology at greater sediment depth 
than the previous investigations and for the first time com- 
pares community structure between a high salinity and low 
salinity site. 

MATERIALS AND METHODS 
SITE DESCRIPTION 

Two sites in Palaeochori Bay [Rocky Point (RP) and Twinkie 
(TW)] were investigated in this study; both are ~100m offshore 
in ~4.5 m deep water (Figures 1C, 2A,B). Preliminary data indi- 
cate nearly identical pH and temperature profiles, although the 
RP site was demonstrably more saline than the TW site. The RP 
site features concentric rings of colored sediments, with orange in 
the center, transitioning to white and then brown (Figure 2A). 
The TW site is located ~8 m due east of the RP site and con- 
sisted primarily of white mats with some gas bubbles at the north 
side of the patch. Samples for this study were collected from the 
white areas from both sites ("P-W" in Figure 2A and "TW" in 
Figure 2B). 

Field sampling protocols 

Pore fluids and sediments for geochemical and microbiologi- 
cal analyses were collected by SCUBA. Sediments were cored 
with polycarbonate tubes (Arthur Kriiger GmbH) and capped 
underwater with rubber end caps. Pore fluids were either col- 
lected in situ into a BD™ (Becton, Dickinson and Company) 
60 mL syringe from a tube inserted into the sediments, or 
extracted through holes in the core tube once onshore using 
10 mL syringes attached to rhizons (long filters which can 
be inserted into the core; Rhizosphere) following Seeberg- 
Elverfeldt et al. (2005). Sediment cores were sliced for subse- 
quent microbiological analyses. Each core slice was aseptically 
placed into sterile bags and frozen on dry ice. Samples were 
then shipped on dry ice back to the laboratory and kept frozen 
at -80°C. 

Samples for free gas geochemistry were obtained by invert- 
ing a serum bottle under water and holding it over the streaming 
gas bubbles. Once the gas had completely filled the bottle, it was 
capped and crimped underwater with little to no seawater inside 
the bottle. Pore fluids for dissolved gases were sampled at 10 cm 
depth in the sediments by inserting a tube into the sediments, and 
filling a syringe. Subsequently, these pore fluids were put into N2- 
filled 60 mL serum bottles until a reverse meniscus formed, then 
capped and crimped. 
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Geochemistry 

Temperatures were measured in situ using a temperature probe 
in a custom-built underwater housing (constructed at the Max 
Planck Institute for Marine Microbiology, Bremen, Germany). 
The pH was measured on shore using a WTW pH meter 3210 
with Mic-D electrode with temperature compensation. Analytical 
uncertainties were approximately ±0.1°C for temperature and 
±0.1 for pH. Pore fluid concentrations of Fe 2+ , NO^, NO^ , 
and NH3 were analyzed in the field using a HACH spectropho- 
tometer. Samples for free (CO2, O2, N2, He, H2, CO, and CH4) 
and dissolved (He, H2, O2, N2, CO, CH4, and CO2) gases 
were shipped on ice to the Istituto Nazionale di Geofisica e 
Vulcanologia, Palermo (Italy), for analysis. Pore fluid samples 
for anion analysis (Br, CI, and SO4 - ) were filtered in the field 
(0.2 |xm), placed on dry ice, and kept frozen until measure- 
ment in the laboratory. Samples for analysis of major seawa- 
ter cations (Na, Ca, K, B, Sr), minor elements (Si, Ba, Mn, 
Fe), and trace elements (Li, Rb, Cs) were preserved in the 
field by filtering (0.2 u,m) and acidification with 0.1% ultrapure 
HNO3. 

LABORATORY PROTOCOLS 

Geochemistry 

Major anions were analyzed using a Dionex ion chromatogra- 
phy system, whereas major cations were measured by induc- 
tively coupled plasma-optical emission spectrometry (ICP-OES; 
Perkin-Elmer Optima 3300). Lithium, rubidium, and cesium 
were measured by inductively coupled plasma-mass spectrome- 
try (ICP-MS; High-resolution double-focusing ICP-MS Thermo 
Finnigan Element 2). Free (CO2, O2, H2, CO, and CH4) and dis- 
solved (H2, O2, CO, CH4, and CO2) gases were measured on a 
Perkin Elmer 8500 gas chromatograph equipped with a double- 
detector (TCD-FID) using argon as carrier gas (Italiano et al., 
2009). 

Microbiology 

DNA extraction, PCR, and Cloning. Bulk environmental 
(genomic) DNA was extracted from ~0.5 g of sediment accord- 
ing to the protocol outlined in the MPBio FastDNA® spin kit, 
and stored at — 20°C. Amplification of 16S rRNA genes for 
both Archaea and Bacteria was performed by polymerase chain 
reaction (PCR) using primer sets 27F/1492R and 21F/1391R, 
respectively (Lane, 1991), on a Hybaid PCR express thermocycler. 
A total reaction volume of 20 or 50 \iL was prepared, containing 
1 x PCR buffer, 0.25 mM each dNTPs, 0.5 |xM each forward and 
reverse primer, 0.5-2 (jlL DNA template depending on sample 
concentration, and 5 U Taq DNA polymerase (5 PRIME). DMSO 
was added to some samples to increase reaction yield. The 
final volume was adjusted with autoclaved milli-Q water. The 
amplification was initiated with denaturation at 95° C for 5 min, 
followed by denaturation (at 95°C for 1 min), annealing (at 
55°C for Archaea, 52°C for Bacteria, for 1 min), and elongation 
(at 72° C for 1 min) for 30-35 cycles. A final extension step was 
carried out at 72°C for 10 min. Agarose gel electrophoresis was 
performed to verify amplification. 

Degenerate primer sets used for arsenite oxidase functional 
gene (azoA-like) amplification were AOX-F-A2 (5'-TGC 



ATCGTCGGCT GYGGNTAY-3') and AOX-R-E2 (5'- 
TTCGGAGTTATAG GCCGGNCKRTTRTG-3') (Zargar et al, 
2012) on a Hybaid PCR express thermocycler. These primers 
target the aoxB gene sequence of the arsenite oxidase operon and 
typically result in ~670 bp product. A total reaction volume of 
50 |xL was prepared, containing 1 x PCR buffer, 0.25 mM each 
dNTPs, 0.5 (jlM each forward and reverse primer, 0.5-2 u,L DNA 
template depending on sample concentration, and 5 U Taq DNA 
polymerase (5 PRIME). DMSO was added for some samples to 
increase reaction yield, and the final volume was adjusted with 
autoclaved milli-Q water. The amplification was initiated with 
denaturation at 95° C for 1 min, followed by denaturation (at 
95°C for 0.5 min), annealing (at 57°C for 0.5 min), and elonga- 
tion (at 72°C for 1 min) for 35 cycles. A final extension step was 
carried out at 72° C for 5 min. Agarose gel electrophoresis was 
performed to verify amplification. 

PCR products were cloned separately using the pCR4-TOPO® 
Plasmid Vector Kit with One Shot TOP10 chemically compe- 
tent Escherichia coli cells. Volumes of 25, 50, and 100 u,L of each 
transformation reaction were spread on selective LB agar plates 
containing ampicillin and incubated overnight at 37°C. White 
colonies were transferred with sterile toothpicks into 2.5 mL, 96 
well microtiter plates containing a final volume of 200 u,L of LB 
medium. After an overnight incubation at 37°C, 100 u,L aliquots 
of liquid cultures were pipetted into a 96 well standard microplate 
with 20 % glycerol and shipped frozen to Beckman Coulter for 
Sanger sequencing. 

16S rRNA gene clone library analysis. Forward and reverse 
sequence fragments were trimmed and assembled using Geneious 
version 5.6.3, created by Biomatters (http://www.geneious.com/). 
Contigs were checked for chimeras using Mothur (Schloss et al., 
2009) and Bellerophon (Huber et al, 2004). Sequences were clas- 
sified using Mothur 's Bayesian classifier and the RDP training set 
version 9 (Wang et al., 2007). Sequences were also aligned, com- 
pared, and clustered into operational taxonomic units (OTUs) 
using the default settings in Mothur. Representative OTUs at the 
97 % sequence identity were selected and used to search for simi- 
lar sequences in NCBI's nr/nt database using nucleotide BLAST 
(http://blast.ncbi.nlm.nih.gov/Blast.cgi). Organisms with com- 
plete genomes, closest cultured representatives, or sequences iso- 
lated from relevant hydrothermal or arsenic-rich environments 
were added as references to the 16S rRNA and awA-like phy- 
logenetic trees. The 97% OTU representatives and the database 
reference sequences where realigned in Mothur, and imported 
into Geneious. The PHYML plugin (Guindon and Gascuel, 2003) 
was used to construct a maximum likelihood tree for bacterial and 
archaeal sequences. 

Several richness and evenness calculations were performed 
using their respective commands in Mothur. These include cal- 
culations for rarefaction curve data, Sobs (the observed richness 
reported, OTU), chao (the Chaol estimator for richness), and 
the simpsoneven and shannoneven (both index-based measures 
of evenness), all at 97% cutoff. 

For az'oA-like gene sequences, assembly and vector trim- 
ming was also done in Geneious. Trimmed sequences were 
searched against Ref_Seq using BLASTx. Sequences that hit 
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arsenite oxidase subjects with an e-value lower than e-40 were 
used to build the tree. Selected bacterial az'oA-like genes where 
taken from Genbank and also added to the tree for reference. 
Sequences were aligned and compared in Geneious using both 
nucleotide and amino acid translation. The aioA tree was built 
using PHYML. 

SUBMISSION OF NUCLEOTIDE SEQUENCE AND ACCESSION NUMBERS 

Sequence data for 16S rRNA bacterial and archaeal clone 
libraries were submitted to GenBank under accession num- 
bers KF278514-KF278557 and KF278485-KF278513, respectively. 
Partial nucleotide sequences for the arsenite oxidase genes were 
submitted under accession numbers KF303544-KF303564. 

RESULTS 
GEOCHEMISTRY 

Table 1 presents the geochemical data, including temperature, 
pH, Fe 2+ , NO3", NO2 , NH3, and major and minor elements, 
collected from RP and TW sites. For comparison, seawater val- 
ues, as well as a depth profile from a control site, were included. 
Temperature and pH were very similar between sites, but Mg 2+ 
and SO4 - concentrations were depleted in both RP and TW sam- 
ples. Furthermore, major cations, such as Na, Ca, K, B, and Sr, 
were enriched in the RP samples, but depleted in the TW sam- 
ples, relative to seawater. Major anions, such as CI and Br, are 
also enriched in RP samples but depleted in TW samples, rela- 
tive to seawater. Most minor elements, including Si, Li, Rb, Fe, 
and Ba, were enriched in both fluids, although only slightly so in 
the TW fluids. Mn was only enriched in the RP samples. Arsenic 
concentrations followed a different pattern compared to all other 
elements. Both sites were highly enriched in this element, but in 
this case, the TW site was much more so, with a maximum con- 
centration of 31.7 (jlM compared to 9.4 jjlM at RP. Depth profiles 
of temperature, pH, and concentrations of CI and As are also 
shown in Figure 3. Note that the values of temperature and pH 
for RP and TW were similar with depth, but values for CI and As 
differed markedly between sites (Figure 3). Table 2 presents the 
dissolved and free gas data collected for this study. Dissolved gases 
were sampled at 10 cm depths at equivalent temperature and pH. 
Both RP and TW free gas was mostly CO2, followed by N2, O2, 
H2, CH4, He, and CO, although H2S concentrations in the free 
gas phase were not analyzed. Dissolved gases followed a similar 
trend. 

MICROBIOLOGY 

Richness and evenness 

The total number of clones generated for Bacteria per site is 
as follows: RP 0-1.5 (51), RP 3-4.5 (41), RP 9-10.5 (53), TW 
0-1.5 (56), TW 3-4.5 (23), and TW 9-10.5 (27). The total num- 
ber of clones generated for Archaea per site was: RP 0-1.5 (27), 
RP 3-4.5 (34), RP 9-10.5 (47), TW 0-1.5 (40), TW 3-4.5 (35), 
and TW 9-10.5 (66). Rarefaction analyses of the 16S rRNA gene 
sequences for Bacteria (Figure Al A) and Archaea (Figure A1B) 
indicated that the RP 0-1.5 sample displayed the highest diver- 
sity, with 22 0TUs out of 51 sequences. The steep slope for the 
RP 0-1.5 sample indicates that a large fraction of the species 
diversity remains to be discovered. The same is true for the RP 



3-4.5 layer, which had 13 0TUs out of 41 sequences. The RP 
9-10.5 sample had the lowest bacterial diversity for all sam- 
ples analyzed. However, the TW 0-1.5 sample had the lowest 
bacterial diversity of all TW samples, with 12 OTUs out of 56 
sequences. The TW 9-10.5 layer diversity was similar to TW 
0-1.5, only slightly more diverse, and a slight development of 
an asymptotic curve. The TW 3-4.5 layer had the highest bac- 
terial diversity of all TW samples, with a slope similar to the 
RP 0-1.5 sample. The RP site showed a clear decrease in bac- 
terial diversity vs. depth, although the TW site did not follow 
this pattern. Hydrothermal environments, compared with many 
other marine environments, are often dominated by relatively few 
taxa (Polz and Cavanaugh, 1995; Schrenk et al., 2004). Therefore, 
in theory, fewer clones may be necessary for adequate sam- 
pling and evaluation of diversity. However, it is important to 
note that these curves may be most useful for assessing relative 
differences in the diversity of major groups between sites and 
depths. 

For the archaeal OTUs, both sites showed increasing diversity 
vs. depth (Figure A1B). Archaeal OTU diversity for the TW 0-1.5 
layer was the lowest, sequences from TW 3-4.5 slightly more so, 
with the TW 9-10.5 layer having the highest diversity. The RP 0- 
1.5, RP 3-4.5, and RP 9-10.5 samples followed a similar trend, 
with increasing archaeal diversity vs. depth. Each of the archaeal 
rarefaction curves level off to the right of the plot, suggesting the 
archaeal diversity may have been sampled adequately. 

Richness indices for Bacteria and Archaea for each site vs. 
depth are presented in Table 3. The Chaol index indicated that 
generally the bacterial communities were more diverse in the 
surface layers of the RP site relative to any other sample, and 
both RP and TW bacterial richness decreased with depth. Both 
the Shannon index and Simpson index indicated that the most 
diverse samples were RP 0-1.5 and TW 3-4.5, and that the RP 
site bacterial richness decreases vs. depth, whereas the TW site has 
the highest bacterial richness in the 3-4.5 depth sample. Each of 
these three indices generally indicated that archaeal richness was 
highest in the deepest samples with overall richness increasing vs. 
depth. 

Bacterial affiliations 

The diversity and distribution of bacterial 16S rRNA clones at the 
phylum (Figure 4A) and class (Figure 4B) levels were investigated 
for three depths at RP and TW. The RP 0-1.5 cm sediments were 
dominated by the phylum Proteobacteria (40 out of 51 clones), 
followed by representatives from the Firmicutes (6), Bacteroidetes 
(2), Actinobacteria (1), Aquificae (1), and Deferribacteres (1) 
(Figure 4A). Within the Proteobacteria, the Epsilon class domi- 
nated (Figure 4B), with nearly all clones identified as Arcobacter 
spp. Minor classes include Deltaproteobacteria (7 clones), Bacilli 
(6), and 1-2 clones each affiliated with Actinobacteria, Aquificae, 
Flavobacteria, and Deferribacteres (Figure 4B). In the 3-4.5 cm 
layer, the Firmicutes phylum dominated (33/41), followed by 
Proteobacteria (6), and one clone each of Bacteroidetes and 
Actinobacteria (Figure 4A). The Firmicutes consisted almost 
entirely of Bacillus spp. (100% confidence); the Proteobacteria 
were composed of the Epsilon (4) and Gamma classes (2) 
(Figure 4B). The deepest layer investigated at RP (9-10.5 cm) 
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FIGURE 3 | Depth profiles of temperature, pH, chloride and arsenic 
concentrations at the RP and TW sites. Note (1 ) similarities between the 
RP and TW site with respect to T and pH, and (2) differences between these 


sites with respect to CI and As concentration. Seawater values are indicated 
with a filled circle, RP data are indicated with a dark gray filled diamond, and 
TW data are indicated with a light gray filled triangle. 


Table 2 | Free and dissolved gas collected for this study*. 


FREE C0 2 (%) 


0 2 (%) 




N 2 (%) 


He (ppm) 


H 2 (ppm) CO (ppm) 


CH4 (ppm] 


RP 92.5 


0.13 




0.67 




7.0 




11450 


0.7 


916 


TW 93.5 


0.39 




1.2 




11.7 




14635 


2.2 


1890.0 


[mL(STP)/L # ] 


Dissolved C0 2 


o 2 




N 2 




He 




H 2 


CO 


CH 4 


RP 50.9 


0.0256 




2.4 




**bdl 


0.0077 0.00029 


0.0058 


TW 245.9 


0.0109 




3.2 




0.00027 


0.0004 0.00027 


0.0303 



'Note: TW free and RP and TW dissolved gases are averages (n = 3; free gas for RP was only collected once). 
**bdl, below detection limit. 

* Dissolved gas data reported as mL of gas at standard temperature and pressure per liter of water. 



was also dominated by the Firmicutes phylum (51/53; Figure 4A), 
with 47 clones identified as Bacilli (Figure 4B), predominantly 
within the genus Bacillus. 

At TW in the 0-1.5 cm layer, Bacteroidetes dominated (37/56; 
Figure 4A), with 36 of those affiliating with the Flavobacteria 
(Figure 4B). The other phyla include Proteobacteria (16), 
Actinobacteria (2), and Thermotogae (1) (Figure4A). Within 
the Proteobacteria, 14 affiliated with the Epsilon class (genus 
Arcobacter), and 2 with the Gamma class (Figure 4B). The 
TW 3-4.5 cm layer was dominated by the Aquifiacae (13 of 
23 clones; Figures 4A,B), especially Thermosulfidibacter spp. 
(often with 100% confidence). The other phyla in this sam- 
ple included Actinobacteria, Planctomycetes, Proteobacteria, 
Thermodesulfobacteria, and Thermotogae (Figure 4A). The 
deepest (9-10.5 cm) sediment sample at TW was dominated by 
the Proteobacteria, followed by Bacteroidetes, Thermotogae, and 
Chloroflexi (Figure 4A). Within the Proteobacteria, the Epsilon 
class (Arcobacter spp.; 18 of 27) dominated, followed by the Beta 
(3) and Gamma classes (2) (Figure 4B). 

The bacterial 16S rRNA phylogenetic tree (Figure 5) includes 
the OTUs from our 16S rRNA clone libraries and some of their 
closest relatives, both isolates and other uncultured clones. Three 



OTUs represented > 5 clones each, which are noted in the Figure 5 
with a green star (173 total clones), and 50 OTUs represented <5 
clones (78 total clones). OTU 23 represents 80 clones from the 
RP site (4 in the top layer, 29 in the middle, and 47 in the deep- 
est layer). This OTU closely affiliates with the Firmicutes class 
{Bacillus spp.; Figure 5). OTU 30 (54 clones) affiliates with the 
Epsilonproteobacteria and was identified in all samples, except RP 
9-10.5 cm. OTU 39 (39 clones) affiliated with Bacteroidetes and 
was identified almost exclusively in the TW 0-1.5 cm layer. 

Archaeal affiliations 

The diversity and distribution of archaeal clones at the three 
depths of RP and TW sediments is shown in Figures 6A 
(phylum) and B (class). At both sites, Euryarchaeota dom- 
inated in the surface layers (Figure 6A). At RP, most of 
the clones affiliated with the Archaeoglobi (16/27; specifically 
the genus Ferroglobus) and Thermococci (8/27; specifically the 
genera Palaeococcus and Thermococcus) (Figure 6B). At TW, 
clones affiliating with Thermoplasmata {Euryarchaeota) dom- 
inated (31 of 40 clones; Figure 6B; specifically the genus 
Thermogymnomonas). Two clones affiliated with Thermoprotei 
{Crenarchaota) . At both sites, with increasing depth, the phyla 
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FIGURE 4 | Diversity and distribution of clones within bacterial phylum 
(A) and class (B) clone libraries for three depths of RP and TW 
sediment layers investigated in this study. Each color represents the 
percentages of each taxon in the total respective clone library, as 
determined by screening of full length 16S rRNA gene sequences. Depths 
are 0-1.5, 3-4.5, and 9-10.5 cm, and indicated at the bottom of site 
designation. 



Table 3 | Microbial diversity richness indices for bacterial and archaeal 
depth profiles from RP and TW sites in Palaeochori Bay. 



Kingdom 


Site 


obs 


Chao 


Shannon 


Simpson 




RP 0-1.5 cm 


22 


73 


2.2 


0.2 


s 


RP 3^1. 5cm 


13 


79 


1.3 


0.5 


E 


HP 9-10. 5cm 


4 


5 


0.5 


0.8 


AC 


TW D-1 R rm 


12 


30 


1 4 


0 4 


ta 














TW 3-4.5 cm 


11 


19 


2.2 


0.1 




TW 9-10.5 cm 


8 


10 


1.4 


0.4 




RP 0-1.5 cm 


9 


10 


1.9 


0.2 


s 


RP 3-4.5 cm 


14 


19 


2.3 


0.1 


1 


RP 9-10.5 cm 


14 


19 


2.2 


0.1 
















TW 0-1.5 cm 


10 


12 


1.3 


0.5 














TW 3-4.5 cm 


11 


26 


2.1 


0.1 




TW 9-10.5 cm 


20 


42 


2.4 


0.1 



'Sobs, number of OTUs identified. 



changed from predominantly Euryarcharota to more and mostly 
Crenarchaeota. The middle layer at RP (3-4.5 cm) was domi- 
nated by Archaeoglobi (Euryarchaeota), followed by Thermoprotei 
(Crenarchaeota; Figure 6B). The Archaeoglobi were again affiliated 
with genus Ferroglobus, while the Thermoprotei genera included 
Staphylothermus, Pyrodictium, and Thermofilum. The deepest RP 
site (9-10.5 cm) was dominated by Thermoprotei (Crenarchaeota; 
Figure 6B), with 42 of 47 clones affiliating with 5 genera includ- 
ing Thermofilum (25) and Staphylothermus (10). In the mid- 
dle TW layer (3-4.5 cm), Euryarchaeota affiliated predominantly 
with Archaeoglobi (15; Ferroglobus sp.), and Thermoplasmata (9; 
Thermogymnomonas sp.), while the Crenarchaeota were entirely 
composed of Thermoprotei (5) (Figure 6B; specifically the gen- 
era Hyperthermus and Thermofilum). The TW 9-10.5 sequences 
affiliated with Thermoprotei (Crenarchaeota; 30 clones; genus 
Thermofilum), Archaeoglobi (Euryarchaeota; 11 clones; genus 
Ferroglobus) and Thermoplasmata (Euryarchaeota; 22 clones; 
genus Thermogymnomonas) (Figure 6B). Note, only one clone in 
all the libraries (RP 9-10.5) affiliated with the Korarchaeota. 

The 16S rRNA archaeal phylogenetic tree (Figure 7) includes 
the OTUs from our 16S rRNA clone libraries and some of 
their closest relatives, other clones, and related isolates. Twelve 
OTUs represented >5 clones each, and are indicated by a 
green star in Figure 7 (212 total clones), and 20 OTUs repre- 
sented <5 clones (37 total clones). Primary OTUs affiliated with 
Thermoprotei (Crenarchaeota), and Thermococci (Euryarchaeota; 
Figure 7). Many sequences were <90% similar to other archaeal 
representatives found in public databases and therefore do not 
cluster with any major group. A considerable fraction of our 
archaeal sequences (~20%) contained large ~26-80bp "inserts," 
possibly introns (Burggraf et al., 1993). The OTU sequences 
containing these inserts (OTUs 8, 1, 30, 5, 3, 27, 2, 33, and 
4) passed comprehensive quality and chimera checks, and are 
therefore included in the phylogenetic tree, and indicated with 
an asterisk, "*" (Figure 7). BLAST hits indicate they are related 
to the Euryarchaeota (Archaeoglobi or Thermoplasmata), and 



Crenarchaeota (Thermoprotaeles or Desulfurococcales). Five of the 
top six OTUs with the highest number of clones contained these 
inserts (8, 1, 30, 3, and 27). OTU 27 (39 clones) was found 
in all samples except RP 0-1.5 (Figure 7), and is identifiable to 
Thermoprotei according to the RDP classification. Three OTUs 
(8, 30, and 1; 38, 38, and 13 clones, respectively) were distantly 
related to Archaeoglobi or Thermoplasmata. OTU 12 (25 clones) 
poorly affiliated with Archaeoglobi (38% ID), which results in it 
being separate from the other groups in the phylogenetic tree 
(Figure 7, top left). 

Arsenic functional genes 

A phylogenetic tree of the aioA-like functional genes (AFGs) 
is shown in Figure 8, and includes sequences from our clone 
library, related isolates and other uncultured clones. Sequences 
affiliated with those detected in Beta (Ralstonia and Burkholderia 
genera), and Alpha classes of the Proteobacteria (Polymorphum, 
Roseovarius and Bradyrhizobium genera) cluster into two distinct 
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FIGURE 5 | Maximum-likelihood phylogenetic tree of closest relatives 
found for representatives from the 16S rRNA bacterial clone libraries. 

OTUs are defined by a 97% similarity cutoff. Symbols next to the OTUs from 
this study represent the samples from which they were recovered, as 



indicated by the sample key at the bottom left. The number of clones 
belonging to each OTU per site is in parentheses. Bootstrap values of >50% 
are shown at nodes. Sulfitobacter mediterraneus was used as an outgroup. 
Green stars indicate the OTUs with the highest number of clones. 
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FIGURE 6 | Diversity and distribution of clones within archaeal phylum 
(A) and class (B) clone libraries for three depths of RP and TW 
sediment layers investigated in this study. Each color represents the 
percentages of each taxon in the total respective clone library, as 
determined by screening of full length 16S rRNA gene sequences. Depths 
are 0-1.5, 3-4.5, and 9-10.5 cm, and indicated at the bottom of site 
designation. 



lineages (Figure 8). Most of the AFGs associated with Beta class 
had no closely related cultured representative. A total of 105 
clones across 22 taxa contained the AFG. Of these, 93 were from 
the TW site, with 57, 24, and 12 clones from the surface, middle, 
and deepest layers, respectively. At the RP site, the AFG was found 
only in the surface (5 clones) and middle (7 clones) layers. 

DISCUSSION 
MICROBIAL DIVERSITY 

Our 16S rRNA gene sequence data indicate that bacterial and 
archaeal communities in the hydrothermally influenced sedi- 
ments of Palaeochori Bay differ significantly when comparing the 
same depths of the RP high salinity and TW low salinity sites, 
although they have approximately the same temperature and pH 
(Figures 4-8, Table 1). Archaeal and bacterial class distribution 



between each site and vs. depth is shown in Table 4, and indi- 
cates that several classes were exclusive to either the RP or TW site. 
For example, comparing the 0-1.5 cm sediment depth, Aquificae, 
Deferribacteres, Bacilli, and Deltaproteobacteria were only present 
in the RP site, whereas Sphingobacteria and Thermotogae were 
only present in the TW sample (Table 4; Figures 4A,B). The same 
is true for the other sediment horizons, with distinct bacterial 
affiliations only occurring at one or the other site. Of the major 
bacterial groups, the Firmicutes were only present in the RP sam- 
ples, while none were present in any of the TW sample depths. 
The 9-10.5 layer has zero overlapping bacterial classes between 
RP and TW sites (Table 4). Betaproteobacteria were identified in 
the TW sample from this depth, the only time this group was seen 
in any of the samples. 

Differences in the archaeal communities in surface samples 
are also quite distinct; e.g., Archaeoglobi is dominant at the 
RP site while Thermoplasmata was dominant at the TW site 
(Figures 6A,B). As was the case for the Bacteria, there exist some 
groups at RP, which were not found at the TW site, and vice 
versa (Table 4). For example, the RP surface sample did not con- 
tain Thermoprotei, a dominant archaeal group for the TW site. 
In the deepest layer, the RP site contained Thermococci, but no 
Thermoplasmata (Table 4). These results suggest that the major 
archaeal and bacterial groups in Palaeochori Bay are distinctly 
different in the high vs. low salinity sites. 

THE ROLE OF PHASE SEPARATION 

Subcritical phase separation is a ubiquitous process in hydrother- 
mal systems, but its effect on microbial communities is rarely 
evaluated. Investigations of the Iheya North Field back-arc 
hydrothermal system south of Japan, suggested that subseafloor 
phase-separation and segregation may influence the supply of 
energy and carbon to vent-associated chemolithoautotrophs and 
related microbial communities (Nakagawa et al., 2005). When 
comparing microbial communities of gas-depleted "normal" 
fluids (i.e., those related to the remaining fluid phase after 
phase separation has taken place) to gas-enriched fluids (i.e., 
those which are related to the vapor phase; German and Von 
Damm, 2003), chemolithotrophs associated with gas-dependent 
energy metabolism, such as hydrogenotrophic methanogenesis, 
were more abundant. Another study investigated differences in 
microbial communities in phase separated hydrothermal fluids 
from the Yonaguni Knoll IV hydrothermal field in the south- 
ern Okinawa Trough (Nunoura and Takai, 2009). While fluids 
showed variability in the physical and chemical compositions 
(e.g., gas content), the microbial communities were relatively 
similar. There, chemolithotrophic hydrogen oxidation was the 
primary metabolism, and it was suggested that similarities in 
microbial communities could have been due to elevated con- 
centrations of hydrogen in both types of fluids. Thus, microbial 
diversity may be influenced by a change in the concentration of 
potential electron donors as a result of phase separation (i.e., dis- 
solved gases such as H2), not necessarily from salinity differences 
alone, although halotolerance cannot be ruled out. As we discuss 
below, metabolisms based on dissolved gases such as H2 are pos- 
sible at each site investigated here, but do not seem to play a 
dominant role in influencing major groups. 
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FIGURE 7 | Maximum-likelihood phylogenetic tree of closest relatives 
found for representatives from the 16S rRNA archaeal clone 
libraries. OTUs are defined by a 97% similarity cutoff. Symbols next to 
the OTUs from this study represent the samples from which they were 
recovered, as indicated by the sample key in Figure 5. The number of 



clones belonging to each OTU per site is in parentheses. Bootstrap 
values of >50% are shown at nodes. Methanomethylovorans hollandica 
was used as an outgroup. Sequences containing inserts are noted with 
an "*." Green stars indicate the OTUs with the highest number of 
clones. 



Our geochemical results indicate that salinity and arsenic 
concentrations are among the most variable parameters when 
comparing RP with TW pore fluid compositions (Tables 1, 2; 
Figure 3). Free gas data indicate that O2, N2, He, H2, CO, 
and CH4 were only slightly higher for TW compared to RP 
(Table 2). Dissolved gas data suggest that no significant dif- 
ferences between sites exists for O2, CO, and H2 concentra- 
tions; the dissolved CH4 is slightly higher in the TW pore flu- 
ids (Table 2). For the redox sensitive species analyzed in this 



study, dissolved Fe 2+ , NO^, and NO2 were slightly higher, 
while NH3 was slightly lower, in fluids from the TW site 
(Table 1). Hydrogen sulfide concentrations can be very elevated 
(up to 3 mM) in high salinity fluids, but were essentially equal 
(~0.25 mM) in both high- and low-salinity hydrothermal flu- 
ids analyzed for this study (Druschel et al., unpublished data; 
Table 1). These geochemical parameters suggest that the phase 
separation process at Milos may create differences in concen- 
trations of dissolved gases, as well as some other important 
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FIGURE 8 | Maximum-likelihood phylogenetic tree of closest relatives 
found for representatives from the aioA-Wke functional gene (AFG) clone 
libraries. AFGs are defined by a 97% similarity cutoff. Symbols next to the 



OTUs from this study represent the samples from which they were recovered, 
as indicated by the sample key in Figure 5. Bootstrap values of >50% are 
shown at nodes. Marinobacter santoriniensis was used as an outgroup. 



redox couples, which could ultimately be contributing to the 
variability in microbial communities when comparing high- 
vs. low-salinity sites. However, we suggest that these differ- 
ences are not significant, and therefore differences in microbial 
diversity may be influenced more by salinity and/or arsenic 
concentrations. 

DOMINANT BACTERIAL METABOLISMS 

The genus Arcobacter, which dominates the RP 0-1.5 cm sample, 
can grow at moderately thermophilic and halophilic conditions, 
utilizing nitrogen or sulfur, including sulfide oxidation (Wirsen, 
2004; Donachie et al, 2005; Pati et al, 2010). For example, A. 
nitrofigilis, isolated from the roots of a salt marsh plant, can fix 
nitrogen (Pati et al., 2010). "Candidatus Arcobacter sulfidicus" 
has been shown to produce large white sulfur mats in hydrother- 
mal environments (Wirsen, 2004), and related organisms maybe 
responsible for the large white mats observed in Palaeochori Bay 
(Sievert et al., 1999). The white sulfur mats contain elemental sul- 
fur, which may be the end product of sulfide oxidation. However, 
a comparison of our representative 16S rRNA gene sequence 



(OTU 30; Figure 5) against this and other Arcobacter species indi- 
cated only a 91-93% sequence similarity (BLAST comparison 
Max ID). 

In the deeper RP sediments, Firmicutes dominate, in particu- 
lar, Bacillus sp. This group is phenotypically highly diverse, and 
includes obligate aerobes, facultative anaerobes, chemotrophs, 
and heterotrophs, halotolerant members, and organisms that 
tolerate a wide range of temperatures (Graumann, 2012). Our 
representative sequence (OTU 23) is most closely related to 
B. licheniformis DSM 13 (Figure 5) which, along with many 
other members of the Bacilli, can reduce nitrate to N2 (Nakano 
and Zuber, 1998). Other investigations suggest that Mn(II) oxi- 
dation (Dick et al., 2006) and arsenate reduction may also 
be mediated by Bacillus spp. (Switzer et al, 1998). Four new 
strains of Bacillus were isolated from shallow-sea hydrother- 
mal vents — two from Panarea and Vulcano Islands (Italy) 
(Maugeri et al., 2002; Spano et al., 2013), and two from the 
sediments in Palaeochori Bay (Milos) (Sievert et al., 2000a; 
Schlesner et al, 2001). The Milos isolates are Halothiobacillus 
kellyi, a mesophilic, obligately autotrophic, strictly aerobic, 
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Table 4 | Bacterial and archaeal class distribution arranged by site and by depth. 



Depth (cm) 


Bacterial class 


Archaeal class 






RP 


TW 


RP 


TW 


U— I.D 


Actinobacteria 


Actinobacteria 


M P 


/ hermoprotei 




Aquificae 


NP 


Archaeoglobi 


Archaeocjiobi 




Flavobacteria 


Flavobacteria 


Methanomicrobia 


M P 
N r 




Deferribacteres 


NP 


Thermococci 


Thermococci 




Bacilli 


NP 


Therm op iasmata 


Thermopiasmata 




NP* 


Sphingobacteria 








Deltaproteobacteria 


NP 








Epsilonproteobacteria 


Epsilonproteobacteria 








Gammaproteobacteria 


Gammaproteobacteria 








NP 


Thermotogae 






TAR 
o— H .0 


Actinobacteria 


Actinobacteria 


1 hermoprotei 


1 hermoprotei 




Sphingobacteria 


NP 


Archaeoglobi 


Archaeoglobi 




Bacilli 


NP 


NP 


M e thanomicrobia 




NP 


Aquificae 


Thermococci 


Thermococci 




NP 


PI a nctomycetacia 


Therm op iasmata 


Thermopiasmata 




NP 


Deltaproteobacteria 








Epsilonproteobacteria 


Epsilonproteobacteria 








Gammaproteobacteria 


NP 








NP 


Thermodesulfobacteria 








NP 


Thermotogae 






y— i u.o 


Actinobacteria 


NP 


1 hermoprotei 


I hermoprotei 




Chloroflexi 


NP 


Archaeoglobi 


Archaeoglobi 




Bacilli 


NP 


NP 


Thermopiasmata 




NP 


Flavobacteria 


Thermococci 


NP 




NP 


Thermomicrobia 


Candidatus Korarchaeum 


NP 




NP 


Betaproteobacteria 








NP 


Epsilonproteobacteria 








NP 


Gammaproteobacteria 








NP 


Thermotogae 







*NR not present. 



chemolithoautotrophic, sulfur-oxidizing bacterium (Sievert et al., 
2000a), and Filobacillus milensis, a halophilic, faculatively 
anaerobic, endospore-forming heterotroph (Schlesner et al., 
2001). In addition, 39 of 80 aerobic mesophiles isolated from 
Milos by Dando et al. (1998) are phenotypically similar to 
Bacillus sp. 

The most abundant bacterial group in the surface sample at 
Twinkie (TW 0-1.5) was Bacteroidetes. However, the RDP classi- 
fication confidence was quite low even at the class level (average 
74%), suggesting many possibilities for representatives from this 
phylum. The representative OTU for these sequences is OTU 
39, which had 35 out of 39 clones in the TW 0-1.5 cm sam- 
ple. Due to the low confidence identification of this group, it is 
not affiliated with any other main group in the bacterial phylo- 
gentic tree (Figure 5, bottom left). BLAST results most closely 
identified with uncultured clones, making identification of pri- 
mary metabolisms difficult. The phylum Bacteroidetes contains 
among the most abundant microbes in coastal marine waters 
(Alonso et al, 2007). Cultured representatives from this phylum 



are typically heterotrophs, non-thermophilic, and use complex 
organic substrates, either by aerobic respiration or fermenta- 
tion (Kirchman, 2002). They can also be found in hydrother- 
mal environments. For example, Kormas et al. (2006) detected 
Bacteroidetes spp. on a white smoker chimney at 9° N, East Pacific 
Rise. Recently, Sylvan et al. (2012) suggested that Bacteroidetes 
replaced "normal" microbial communities once smoker chim- 
neys become inactive. Sievert et al. (2000b) and Nitzsche (2010) 
both identified Bacteroidetes as dominant members in the white 
areas from Milos sediments, but temperatures were in both cases 
quite lower than for our study (~20°C less). This suggests the 
possibility that Bacteroidetes activity is limited to lower temper- 
ature environments. Temperatures were comparable between the 
RP and TW sites (i.e., <~ 40°C), but salinity was quite differ- 
ent, suggesting the lack of Arcobacter in place of Bacteroidetes at 
the TW site may also be influenced by halotolerance. However, 
as discussed, dissolved gases and other redox couple concentra- 
tions may play a role, as they are slightly elevated in the TW pore 
fluids. 
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The TW 3-4.5 site was dominated by Thermosulfidibacter 
sp., whose closest cultured representative is Thermosulfidibacter 
takaii, a thermophilic, sulfur-reducing chemolithoautotroph 
from a deep-sea hydrothermal vent at the Yonaguni Knoll IV, 
Southern Okinawa Trouph (Nunoura et al., 2008). Growth was 
observed within the temperature range of 55-78°C, and 0.5- 
4.5% NaCl. Temperatures at this sediment depth should be near 
50° C, and salinity should be slightly lower than seawater (Table 1; 
Figure 3). 

The TW 9-10.5 was dominated by Arcobacter spp. 
{Epsilonproteobacteria; 18 out of 27 clones). As discussed, 
the group could be conducting nitrate reduction or sulfide 
oxidation. Comparing the OTU representative sequence (OTU 
5) from TW 9-10.5 to the 16S rRNA gene sequence from A. 
nitrofigilus and "Ca. A. sulfidicus" indicates BLAST max ID of 
93% and 93%, respectively. Indeed, the representative OTUs 
from RP 0-1.5 (OTU 14) and TW 9-10.5 (OTU 5) were only 
96% max ID when compared to each other. This suggests that 
the Arcobacter species in both the RP 0-15 and TW 9-10.5 
samples are only distantly related, and it is probable that they 
are different species or genera with different environmental 
tolerances. 

These data indicate that bacterial communities are strati- 
fied by depth, and that sulfur- and nitrogen-based metabolisms 
may dominate at the RP site, whereas phototrophy, sul- 
fur, and or nitrogen-based metabolisms may be dominant at 
the TW site. Temperature may be the dominant parameter 
which influences microbial community structure as a func- 
tion of depth at both RP and TW sites, as vertical gradi- 
ents are more pronounced for this parameter compared to 
any other (Table 1; Figure 4). Sunlight may also play a role, 
particularly for the surface TW site where Bacteroidetes is 
present. 

DOMINANT ARCHAEAL METABOLISMS 

The dominant archaeal genera for RP 0-1.5 were Ferroglobus 
(Archaeoglobi) and Palaeococcus/Thermococcus (Thermococci). 
These groups become replaced mostly by the Thermoprotei 
(Staphylothermus and Thermophilum) at depth. However, poor 
RDP classification confidence makes identification beyond 
the phylum level difficult. Further complicating matters, the 
OTU chosen to represent the sequences associated with the 
Archaeoglobi (OTU 8) contained a 44 bp insert. BLAST analysis of 
OTU 8 returned no cultured representatives, and identified most 
closely with uncultured clones from the Logatchev hydrothermal 
field (LHF) and terrestrial hydrothermal vents on Ambitle Island, 
PNG. Removing the insert from OTU 8 and conducting a BLAST 
analysis for the remaining sequence indicates Aciduliprofundum 
spp. as the closest cultured match (100% query coverage; 86% 
max ID). This organism belongs to the "Deep-sea Hydrothermal 
Vent Euryarchaeota 2" (DHVE2) lineage (Takai and Horikoshi, 
1999), which are widespread at deep-sea hydrothermal vents 
(Reysenbach et al., 2006; Flores et al., 2012). Unfortunately, rel- 
atively little is known about their distribution and phylogenetic 
diversity (Flores et al., 2012). The only described species from 
the Aciduliprofundum genus is A. boonei T469 T , which is an 
obligate thermoacidophilic sulfur- or iron-reducing heterotroph 
capable of growing from pH 3.3-5.8 and between 55 and 75°C 



(Reysenbach et al., 2006). This organism is included in the 
phylogenetic tree, but does not affiliate closely with our sequences 
(Figure 7). 

The second most abundant archaeal group in the RP 0-1.5 
sample was Thermococci, consisting of both Thermococcus spp. 
and Palaeococcus spp., with ~100% RDP confidence. Our OTUs 
for this group (31, 6, and 11) did not contain insertions, and 
identified most closely with Thermococcus spp. (Figure 7). There 
are numerous species related to this genus; most are consid- 
ered extreme thermophiles, typically identified as organotrophic 
anaerobes which grow above 70° C, sometimes requiring elemen- 
tal S, producing H2S and CO2, and have the ability to use H + 
(Arab et al, 2000). 

The RP 9-10.5 sample was dominated by the Crenarchaeota- 
Thermoprotei, was composed of several genera, which were dom- 
inated by Thermofilum spp., although mostly with very poor 
confidence, followed by Staphylothermus spp., often with 100% 
confidence, and several others. Our OTUs chosen to represent the 
Thermofilum spp. identified by the RDP classification (21, 26, and 
27) most closely affiliated with Thermoprotei (Figure 7). However, 
OTU 26 is within a slightly different cluster, and had much 
higher RDP % confidence than the other two with Thermofilum 
spp. BLAST analysis of this OTU most closely matched with 
Thermofilum pendens, a thermophilic, anaerobic sulfur-respirer. 
OTU 27 contained an insert, and BLAST analysis indicated 
the closest representative was the uncultured clone VulcPlw.113, 
from the shallow-sea hydrothermal vents off Vulcano, Italy. The 
sequence from this clone, as noted, also contains an insert. 
Removing the insert from OTU 27 and conducting a BLAST 
analysis indicates it does not affiliate with any cultured represen- 
tatives, and is most closely related to various uncultured Archaea 
from the Mariana Trough, East Pacific Rise, and the shallow-sea 
vents off Ambitle Island, Papua New Guinea (Meyer-Dombard 
etal., 2011). 

The TW 0-1.5 was by far dominated by the Thermoplasmata 
class, Thermogymnomonas genus. This genus has one het- 
erotrophic species, T. acidicola, and can grow at temperatures in 
the range 38-68°C (optimally at 60°C), at pH 1.8-4.0 (optimally 
at around pH 3.0), and is obligately aerobic and heterotrophic, 
requiring yeast extract for growth (Itoh et al, 2007). These clones 
are represented by OTU 30, which according to RDP, poorly 
affiliated with Thermoprotei. This OTU contained an insert, 
and BLAST analysis indicated affiliation to uncultured clones 
VulcPlw.76, and those from terrestrial vents from Papua New 
Guinea. The closest cultured representative was Aciduliprofundum 
sp. MAR08-237A. BLAST analysis of this sequence with the insert 
removed indicates similar results. 

In the deeper 9-10.5 sample for TW, dominant members tran- 
sition to Thermoprotei (mostly Thermofilum spp.; OTUs 2 and 
3), although abundant Thermoplasmata (Thermogymnomonas 
spp.; OTUs 1 and 5) were also present. Metabolisms for these 
organisms have been discussed previously. Each of these OTUs 
contained inserts. OTUs 2 and 3 identified most closely with 
VulcPLw.76 uncultured clone and previously obtained clone 
sequences from Milos, with no affiliation to a cultured representa- 
tive. OTUs 1 and 5 identified with the isolate A. boonei. OTU 4 in 
this sample identified with Stetteria sp., from the Thermoprotei- 
Desulfurococcaceae family, by the RDP classification. BLAST 
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analysis of this sequenced indicates affiliation with the previously 
described uncultured clones from the shallow-sea vents off 
Vulcano and Milos. 

These results indicate that archaeal metabolisms in both sur- 
face and deeper sediments are dominated by sulfur oxidizers, can 
be heterotrophic, and may also have the ability to conduct Fe oxi- 
dation. These poorly classified groups, particularly those which 
contain an insert, could represent new genera or new higher order 
taxa. 

Although we cannot say for sure these extra nucleotide 
sequence insertions are introns, we believe the extra nucleotides 
in our sequences are real for three main reasons: (1) The quality 
scores for most of these sequences were very high. (2) Multiple 
chimera checks suggested they are not chimeras (see Materials 
and Methods section), (3) They are found not only in our dataset, 
but also in several other datasets from Milos [Brehmer, 2000 
(unpublished data, NCBI database); Nitzsche, 2010], and the 
near-shore hydrothermal vent system of Vulcano (Italy) (Rogers 
and Amend, 2005). The other datasets were obtained from com- 
pletely different laboratories using different methods, making it 
highly unlikely that there is an error. All of the other introns 
found in this study, when blasted independently after exons were 
removed, consistently returned affiliations with clone sequences 
from these previously reported groups as the closest represen- 
tatives. An intron is an additional piece of nucleotide sequence 
found in the genes of most organisms and many viruses, and can 
be located in a wide range of genes, including those that gen- 
erate proteins, rRNA, and tRNA (Itoh et al, 2003). Introns in 
the 16S rRNA gene of Archaea are rare, and primarily restricted 
to the Crenarchaeota, within the families Thermoproteaceae and 
Desulfurococcales (Morinaga et al., 2002; Itoh et al, 2003). 
Although outside the current scope of this paper, it is important 
to note that if the extra nucleotide sequence contained within our 
archaeal sequences for OTU 8, as well as OTUs 1, 30, and 5, are in 
fact introns, it is a significant finding as introns contained within 
the 16S rRNA gene sequences of Euryarchaeota have not, to our 
knowledge, been reported. 

ARSENIC FUNCTIONAL GENE SURVEY 

Elevated arsenic concentrations may clearly influence micro- 
bial community structure and function, since the concentrations 
encountered can be toxic and may also provide additional redox 
couples for microbial metabolism (e.g., arsenotrophy). For exam- 
ple, in reducing hydrothermal fluids, arsenic exists most com- 
monly as arsenite, with a +3 oxidation state (AsIII, F^AsO^). 
However, in surface seawater and other oxidizing environments, 
arsenic is dominantly in the +5 oxidation state (arsenate, AsV, 
HAsO^). Thus, microorganisms can play an important role in 
the arsenic cycle as they catalyze the transformation of differ- 
ent arsenic redox species (Handley et al., 2009; Akerman et al., 
2011; Meyer-Dombard et al., 2011, 2013). The three main types of 
arsenotrophs thus far known to occur are heterotrophic arsenite 
oxidizers (HAOs), chemotrophic arsenite oxidizers (CAOs), and 
dissimilatory arsenate-reducing prokaryotes (DARPs; Oremland 
and Stolz, 2003). Genes involved in these arsenic transformations 
include detoxification via the ArsC system, and energy gain via 
the arsenate respiratory reductase (two subunits, ArrA and ArrB) 
or the arsenite oxidase (two subunits, AioA and AioB, formerly 



referred to as AroA and/or AoxB; Oremland and Stolz, 2003; 
Zargar et al., 2010). Recently, an alternative AsIII-oxidation gene 
was discovered (arxA-type; Zargar et al., 2010). To date, most 
studies of As microbial cycling have focused on various naturally 
or anthropogenically impacted environments, primarily located 
in terrestrial (non-marine) settings. Next to nothing is known 
about As cycling in coastal marine environments. 

At the two sites investigated here, analysis of the AFGs revealed 
affiliation with only two classes within the Proteobacteria phylum, 
although we must keep in mind that lateral gene transfer may 
have occurred (Figure 8). The Alphaproteobacteria containing 
the arsenite oxidase gene were predominantly affiliated with the 
genera Polymorphum, Bradyrhizobium and Roseovarious, while 
the Betaproteobacteria were affiliated with the Burkholderia or 
Ralstonia genera. The primers used for screening of aioA-\ike 
genes were based on a range of groups, and therefore should 
be broadly applicable (i.e., should not be specific to one or two 
groups). For example, primer design relied on the sequence align- 
ment of 16 arsenite oxidase amino acid sequences consisting of 
5 Alpha-, 4 Beta-, 3 Gammaproteobacteria, 2 Deinococci, and one 
Chloroflexi (Zargar et al.,2012). Thus, there could be a bias toward 
Alpha- and Betaproteobacteria, but many published arsenite oxi- 
dizers are in these groups (e.g., the most updated phylogenetic 
tree for arsenotrophs has many of the azoA-like gene contain- 
ing microbes as Alpha- and Betaproteobacteria). There are some 
newer isolates being reported within the Gammaproteobacteria, 
and various genome sequencing projects are also starting to show 
the presence of cHoA-like genes (Handley et al., 2009; Cavalca 
et al., 2013). The primers used here were the best available at the 
time of their design. 

AFG distribution was primarily associated with the TW sam- 
ples, which is somewhat expected given the higher concentra- 
tions of arsenic at this site (Figure 8). Phylogenetic analysis 
revealed that both the Beta- and Alphaproteobacteria relationships 
can be further subdivided into two distinct lineages (Figure 8). 
While some of the Betaproteobacteria were most closely related 
to Ralstonia and Burkholderia genera (AFGs 3 and 4), nearly 
all of the other AFGs associated with Betaproteobacteria had 
no closely related cultured representative (Figure 8). Four main 
clusters exist for these Betaproteobacteria. They consist of (1) 
AFG 6, (2) AFGs 2, 5, 15, 7, 12, 20, and 21, (3) AFG 19, 
and (4) AFGs 8 and 13. BLAST analysis of these AFGs indi- 
cate they are associated with either Beta- or Alphaproteobacteria, 
mostly Burkholderiales, but cannot be constrained beyond the 
Proteobacteria phylum. The Alphaproteobacteria were most closely 
related either to Polymorphum, Roseovarius, or Bradyrhizobium 
spp. AFG 11 was closest to Polymorphum gilvum in the phy- 
logenetic analysis (Figure 8). Only Burkholderiales occurred in 
the low salinity 9-10.5 cm depth in the 16S rRNA gene clone 
libraries. BLAST analysis of these AFG sequences indicates sim- 
ilarity to aioA-Wkt genes from various terrestrial arsenic-rich 
environments. 

These species are only distantly related, and therefore discus- 
sion of their arsenic oxidation strategy is limited. For example, 
Polymorphum gilvum was isolated from a crude oil-polluted saline 
soil in Shengli Oilfield, China, and was able to use the crude oil 
as the sole carbon source, and is not known to conduct arsenic 
metabolism, although its genome contains genes for both arsenite 
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oxidation and arsenate reduction (detoxification) (Li et al., 2011). 
All of the related groups in the AFG analysis can be confirmed 
to contain an arsenite oxidase, but none of them have been 
shown to oxidize or reduce arsenic for detoxification or for their 
metabolism. 

Overlap in AFG and 16S rRNA gene sequence data occurred 
only in the low salinity TW 9-10.5 sediments. For example, two 
OTUs (2 and 3) from the bacterial 16S rRNA data set were 
identified most closely with the Betaproteobacteria, specifically 
Ralstonia sp. and Burkholderia sp. This suggests arsenotrophy 
may be an important metabolism in deep, low-salinity, As-rich 
environments similar to those found in the TW 9-10.5 layer. 

Other molecular investigations of the genes responsible for 
As redox transformations have revealed intriguing evidence for 
As microbial transformations in marine hydrothermal environ- 
ments. For example, Nitzsche (2010) used arsenic primers for 
az'oA-like gene amplification of 0-6 cm pooled sediments from 
Milos, and obtained similar results to our data (i.e., Alpha- 
arid Betaproteobacteria). At a marine shallow- water hydrother- 
mal vent off Ambitle Island, Papua New Guinea, microorganisms 
may be key in the oxidation of As(III) in areas immediately 
surrounding the vent source (Akerman, 2009; Meyer-Dombard 
et al, 2011, 2013). Using PCR primers designed to amplify 
arsenic functional genes, biofilm communities at the site were 
shown to possess the genetic capacity for the oxidation of As(III). 
The aroA-like (now known as azoA-like) clone sequences ana- 
lyzed belonged to relatives of Alpha- and Betaproteobacteria, 
Thermus, and Pyrobaculum, but were only 71-80% similar to 
these sequences available in GenBank, suggesting that many of the 
Ambitle Island As-functional genes are unique. Thus, in marine 
environments, arsenite oxidation may be limited to the Alpha- 
arid Betaproteobacteria, but much more research is necessary, 
particularly cultivation and molecular-based (next generation 
metagenomic) approaches focusing on arsenotrophy. 

CONCLUSIONS AND OUTLOOK 

Our results indicate that differences in archaeal and bacte- 
rial community structure exist at the site as a function of 
not only depth, but also when comparing high vs. low salin- 
ity environments with similar redox potential. A wide range 
of metabolisms seem possible, with affiliated members likely 
involved in the nitrogen or sulfur cycle. However, it should be 
noted that rare members of the community can also be very 
active and function as "keystone species," yet not be detected 



in clone libraries. Differences in bacterial and archaeal commu- 
nities between the RP high-salinity and TW low-salinity sites 
suggest geochemical signatures resulting from phase separation 
may be a controlling factor in hydrothermal vent microbial 
ecology (e.g., salinity, arsenic). Our screening of AFGs respon- 
sible for arsenite oxidation indicate that both the Alpha- and 
Betaproteobacteria may be responsible for arsenite oxidation in 
sediments from both sites, although the most abundant clus- 
ters are from the higher-arsenic TW site. These results are 
similar to previously reported AFG screening in marine envi- 
ronments, suggesting that Alpha- and Betaproteobacteria are 
the key players in marine hydrothermal arsenic metabolism. 
Since these groups are only distantly related to known arsen- 
ite oxidizers, novel genera may be present at the site; much 
more research is necessary including targeted cultivation of 
arsenotrophs. These results also provide incentive for metage- 
nomic investigation of these sites to identify potential novel 
shallow vent species, and confirm arsenite metabolism is present 
primarily in the Alpha- and Betaproteobacteria. Affiliations from 
the arsenic gene screening were often not reflected in the 16S 
rRNA results, suggesting that arsenite may only be an impor- 
tant electron donor in the deeper, low salinity, high arsenic 
sediments. However, this could also be the result of a lack of 
sufficient database representation, relatively shallow sequencing 
depth or potential bias in the AFG primers. Finally, many of the 
archaeal 16S rRNA sequences contained putative introns, includ- 
ing members of Euryarchaeota, which has not been reported 
previously. 
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FIGURE Al | Rarefaction curves showing 16S rRNA OTU richness for each of the sampling sites and depths for this study. Curves were calculated 
using the furthest neighbor approach with a 97% sequence similarity cutoff. Data shown are for Bacteria (A) and Archaea (B). 
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